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The  P.  I.  and  Miss  Decker  st  lied  the  effects  of  introducing  finite 
conductivity  and  compressibility,  as  well  as  allowing  for  a  variety  of 
equilibrium  profiles,  on  the  development  of  the  Kruskal-Schwarzchild 
instability.  The  inclusion  of  finite  conductivity  effects  produced  a  fourth 
order  eigenvalue  equation  for  the  modes.  A  nunerical  solution  of  the  equation 
was  done  and  a  dispersion  relation  obtained.  Miss  Decker  wrote  up  her  work  as 
her  doctoral  dissertation  and  was  granted  the  Ph.  D.  in  the  Spring  of  1986. 
A  draft  of  a  paper  we  intend  to  submit  is  attached. 

The  work  with  Mr.  Boynton  centered  on  a  nunerical  solution  of  the  MHD 
equations  for  a  model  plasma.  A  two  dimensional  FCT  code  was  written  and 
tested  on  a  variety  of  simple  problems  where  the  answers  could  be  obtained  by 
other  means.  A  good  deal  of  effort  was  also  expended  on  developing  graphical 
methods  that  would  allow  good  display  of  the  results.  A  sample  output  figure 
is  attached.  It  is  a  plot  of  pressure  vs.  position  for  a  state  that  evolved 
from  an  initial  pressure  spike.  No  magnetic  effects  are  included  in  the  case 
shown. 

Other  work  was  done.  For  example  the  P.  I.  and  Dr.  Nearing  improved 
the  results  for  the  distribution  of  rail  current  to  the  case  of  arbitrary  time 
dependence  of  the  current  and  acceleration,  as  well  as  a  finite  thickness  for 


the  armature. 


FINITE  CONDUCTIVITY  AND  THE  RAYLEIGH-TAYLOR  INTERCHANGE  INSTABILITY 
Ann  M.  Decker  and  Mannel  A.  Huerta 
University  of  Miami  Physics  Department 


The  literature  of  the  physics  of  fluids  abundantly  documents  the  phenomenon 
that  has  come  to  be  known  as  the  Rayleigh-Taylor1'*,  or  interchange, 

instability.  A  fluid  supported  against  gravity  by  a  lighter  fluid  is  found  to 
be  unstable,  the  potential  energy  of  the  system  being  lowered  by  the  exchange  of 
the  two.  When  the  gravitational  field  is  opposite  the  density  gradient  of  a 
fluid,  any  initial  perturbation  grows  exponentially  in  time. 

The  hydromagnetic  analogue  of  this  instability  was  first  explored  by 

Kruskal  and  Schwartzschild. ’  A  plasma  supported  against  gravity  by  a  magnetic 
field  tangential  to  its  plasma-vacuum  interface  is  subject  to  an  instability  of 
mathematically  identical  form.  A  Rayleigh-Taylor  instability  is  thus  expected 
in  the  case  of  a  plasma  slab  accelerated  by  a  magnetic  field. 

Many  references  exists  describing  the  dispersion  relations  of  special  cases 

of  Rayleigh-Taylor  (R-T),  or  Krnskal-Schwartzschild  (K-S) ,  instabilities  in 

plasmas.  In  the  case  of  a  plasma  slab  described  by  the  laws  of  ideal 

magnetohydrodynamics,  a  surface  deformation  of  wave  vector  k  perpendicular  to 

the  magnetic  field  grows  with  dispersion  u  =  (kg)1^*,  where  g  is  the 

acceleration  of  the  slab.4'*  Tsai,  et  al.*  have  detailed  the  dispersion 

relations  numerically  for  a  compressible  plasma  of  infinite  conductivity. 

Huerta  and  Decker7  have  shown  that  the  growth  of  R-T  instabilities  is  unaffected 

by  the  inclusion  of  finite  conductivity  to  first  order  for  the  case  of  an 

incompressible  plasma.  For  parameters  typical  of  plasma  armatures  in 

rail-launch  devices,  even  the  longest  wavelength  disturbances  would  have  ample 

time  to  develop  fully  at  this  rate.  Such  an  instability  in  the  plasma  arc  would 

10 

destroy  the  type  of  one-dimensional  equilibria  used  in  earlier  arc  models. 

A  fully  developed  R-T  instability  will  result  in  the  breakoff  and  expansion 
of  plasma  packets  away  from  the  main  arc.  Trailing  masses  of  plasma  may  lead  to 
the  formation  of  secondary  arcs  and  turbulent  flow.  Such  effects  have  been 
reported  by  Parker,  et  al.11,  and  are  important  in  that  they  could  limit  the 


efficiency  of  plasma  driven  mass  accelerators.  A  detailed  study  of  the  R-T 
instability  in  a  plasma  slab,  including  viscosity,  compressibility,  finite 
conductivity,  etc.  is  warranted  to  see  what  factors  can  impede  its  growth. 

The  present  discussion  will  include  the  effects  of  finite  conductivity  in 
both  the  equilibrium  and  perturbation  equations.  A  general  mode  equation  will 
be  derived,  which  for  an  infinitely  conducting,  isothermal  plasma  will  reduce  to 
Tsai's  result.  An  equilibrium  profile  consistent  with  a  rail  geometry  and  the 
assumption  of  a  one-dimensional  steady  state  will  be  inserted  into  a 
dimensionless  form  of  the  general  mode  equation.  In  the  limits  of  both  very 
high  and  very  low  conductivity,  the  resulting  equation  depends  only  on  the  ratio 
l g0/K.  where  f  is  the  plasma  tail  length,  g0  is  the  arc  acceleration,  and  K  is 
the  square  of  the  ion  acoustic  velocity. 


The  governing  equations  for  the  plasma  system  will  be  Maxwell's  equations, 
the  momentum  equation,  the  equation  of  continuity,  and  an  equation  of  state. 
Displacement  current  is  neglected.  An  adiabatic  equation  is  assumed, 
p/pr  =  const,  where  as  yet  y  is  unspecified.  The  effects  of  finite,  constant 
conductivity  come  from  E  +  v  X  B/c  =  J/o.  The  equilibrium  state  of  a  neutral, 
stationary  plasma  is  then 


-  J  X  B  + 

c 


V  X  B.  =  — J 


A  plasma  that  has  experienced  a  small  perturbation  from  its  equilibrium 
can  be  described  by  the  variables  B^,  Jj(  E^,  px,  p  ,  taken  to  be  small 
deviations  from  the  equilibrium  values,  B  ,  J  ,  E  ,  p  ,  p  .  These,  and  the 
velocity  v,  are  substituted  into  Maxwell's  equations,  the  force,  and  the 
continuity  equations.  Neglecting  terms  of  second  order  in  the  perturbing 
variables. 


(1) 


Eliminate  various  dependent  variables  in  these  equations  and  the  result  is 
a  general  first  order  perturbation  equation  for  a  plasma  of  arbitrary  initial 
current,  mass  density,  magnetic  field,  acceleration,  and  constant  conductivity. 

a*v1 

P0 — ^  “  V  [v4  •  Vp0  +  YP0V*vj] 

at 

2 

+  5^W7  x  <v®.>]  -  jib T  x  Iv  x  <T  x  «,>]] x  e. 

+  Z.  x  E*  x  <1.  X  E.>  -  s:7  X  (T  X  S,)]  -  (T-P.  I,)*.  (2) 


The  general  equation  above  will  now  be  applied  to  the  specific  cartesian 
geometry  of  interest.  The  equilibrium  parameters  are  as  Bq  =  BQ(y)z, 
£0  c  -J0(y)i,  ^  =  ~g0y,  p0(y) ,  and  p#(y)»  while  the  velocity  and  the 

■*  j  I-* 

perturbation  quantities  will  be  in  the  form  Ax(z,y)  e 

For  this  geometry,  equation  (2)  can  readily  be  analyzed  in  terms  of  its  x 
and  y  components. 


B  a 

— <o*p  y  «=  — p  g  — — v  +  “ — T  (yp  +  )  (V*v  )  —  ■  ■■  -c  ■■■'  B  V*B  1  (3) 
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0)  p  V  *=  O  K  . 

V°  xy  vo  *o  fly  iy 

B*  * 
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16n  a 

The  parenthetical  quantity  can  be  removed  between  these  two  equations  by  taking 
the  derivative  with  respect  to  the  y  coordinate  of  equation  (3),  and  the 
derivative  with  respect  to  x  of  (4)  and  subtracting.  The  result  is 

3vxx  V1X  3pe  r,  *•  1  3v*y  go  a  , 

3y  p  3y  L  *  p  3y  J  3x  a  3x  1 


3y  p#  3y 


a  p  3y 

(I)  r  o  * 


It  should  immediately  be  noted  that  all  terms  involving  a  a  dependence  have 
been  removed,  and  this  result  is  valid  for  a  plasma  of  arbitrary  conductivity. 


To  this  point,  the  results  have  been  applicable  to  a  wide  range  of  plasmas 
obeying  a  very  general  adiabatic  equation  of  state.  Now  the  discussion  will 
narrow  slightly  to  the  consideration  of  plasmas  obeying  the  isothermal 


relations,  p  ■=  Kp,  y  =  1,  as  need  in  other  references  examining  the  same 

4.11 

phenomenon.  Here  K  is  the  square  of  the  ion  acoustic  velocity, 

K  =  (i+lll'T/m^,  where  z  is  the  (constant)  ionization  state  of  the  plasma  ions, 
k'  is  Boltzmann's  constant,  T  is  the  temperature  in  Kelvins,  and  m^  is  the  mass 
of  the  ionic  (or  neutral)  species, 

Powell  has  calculated  the  equilibrium  parameters  under  the  assumption  of 
a  finite  conductivity,  isothermal  plasma.  It  is  found  that  the  equilibrium  is 
described  by  the  relations 

Bp (y)  =  ^pT#(f-  y)  »  P„  (y)  “  ~  “  [f  +  7  -  (f+  X)e  y/J^sqbrk  (6) 

c 

where  Jg  is  constant  and  X  is  the  ratio  of  K,  the  square  of  the  ion  acoustic 
velocity,  to  the  acceleration  go< 

The  governing  relations  for  the  plasma  will  be  scaled  in  terms  of  the 
length  t  of  the  plasma  tail,  and  time  (f/gp)1^*  for  the  projectile  to 
accelerate  this  length.  The  following  transformations  will  be  used: 


a 


(7) 


The  quantities  u,  q,  w,  and  s  are  the  dimensionless  coordinate,  wave  vector, 
perturbation  frequency,  and  conductivity,  respectively.  Furthermore,  the 
magnetic  field  and  mass  density  can  be  expressed  in  terms  of  their 
dimensionless  analogues,  b  and  p. 


arbitrary  initial  current,  density,  and  magnetic  induction  profiles.  We  use 
the  isothermal  equation  of  state  and  express  the  mode  equation  terms  of  R: 


S  =  p0v 


iy 


(9) 
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All  primes  indicate  derivatives  with  respect  to  the  dimensionless  variable  n. 
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(13) 


(14) 


(15) 


Equations  (10)  through  (15)  show  that  the  dispersion  of  q  vs.  w,  to  be 
derived  through  the  integration  of  (10),  can  only  vary  with  s  and  x.  In  both 
the  limits  of  very  high  and  very  low  conductivity,  the  factor  s  is  eliminated, 
so  that  the  behavior  of  the  dispersion  relation  is  characterized  by  the 


l  a 

dimensionless  parameter  x.  Powell  has  already  fonnd  this  to  be  the  case  for 
the  limit  of  large  conductivity,  and  has  done  the  numerical  integration  of  for 
that  limit. 

In  the  limit  of  infinite  conductivity,  the  mode  equation  reduces  to  the 
exact  form  of  the  mode  equation  derived  by  Tsai*  under  the  initial  assumptions 
of  large  conductivity  and  an  isothermal  equation  of  state. 

The  boundary  conditions  at  the  plasma-vacuum  interface  are  largely 
determined  by  two  factors:  the  continuity  of  components  of  the  electromagnetic 
fields,  and  the  fluid  properties  and  physical  constraints  on  the  velocity, 
pressure,  mass  density,  etc.,  and  their  derivatives.  These  constraints  will 
place  certain  restrictions  on  the  derivatives  of  the  velocity,  and  on  the 
derivatives  of  the  function  E  for  this  density  profile.  The  results  of  this 
analysis  are,  in  their  the  dimensionless  form,  and  for  the  specific  density  and 
magnetic  induction  profiles  (8): 


R(  0)  =  0  , 

R'(0)  *»  unspecified  , 

r  *  _  *, 

R”  (0)  =  | 

[a-  -  iws  -  1  -  X  ] 

w 

|R'(0) 

R"(0)  =  [iws 

[2  +  x  +  p.  -  A-]  +  a_[ 

2w  w 

2w*  -  x  - 

l]  +  x[l  +  x  -  w‘] 

|]r'(o) 

The  value  of  R'(0)  will  be  left  as  a  free  parameter,  so  that  the  solution  of 
the  integration  of  equation  (10)  will  be  correct  to  within  a  constant 
multiplicative  factor. 

Equation  (10)  was  integrated  subject  to  the  boundary  conditions  derived 
above  and  the  additional  constraint  that  the  fluid  velocity,  and  hence  R,  must 
vanish  at  the  plasma-projectile  boundary.  Thus,  for  a  given  wave  vector  q,  the 
frequency  can  be  deduced  by  finding  those  values  of  w  that  ensure  R  =  0  at  u  = 
1.  A  Runge-Kutta  algorithm  was  used  for  the  integration. 

As  with  equation  (10),  the  boundary  conditions  vary  only  with  the 
dimensionless  parameters  x  and  s.  The  dispersion  relation  q  vs.  w  have  been 
found  for  several  values  of  these  parameters. 

It  is  of  interest  to  compare  the  relative  sizes  of  s  and  x  for  a  tangible 
model.  Powell11  has  shown  that,  for  small  x,  x  430/2,  where  0  is  the  arc  mass 
to  projectile  mass  ratio.  McNab1*  has  approximated  the  mass  of  the  plasma  arc 


in  the  Rashleigh-Marshall  launch  experiment  to  be  0.1  grams.  For  the  3  gram 
mass  used  in  that  experiment,  the  approximate  value  of  x  is  0.05.  Spitzer** 
gives  the  conductivity  of  a  fully  ionized  plasma  in  terms  of  the  average 
ionization  state,  the  temperature,  and  the  electron  number  density.  For  a  rail 
armature  experiencing  the  peak  current  per  rail  height  5.76  x  1014  statamp/cm 

reported  by  Rashleigh  and  Marshall,  numerical  calculations  by  Powell  and 

io  i  a 

Batteh  give  these  parameters.  Combining  them  with  the  relation  for  the 

conductivity,  the  acceleration,  and  the  definition  of  s,  one  arrives  at  a  value 

of  s  for  the  Rashleigh-Marshall  plasma  armature:  x  =  0.05,  s  =  11.2 

The  relative  values  of  s  and  x  derived  for  the  Rashleigh-Marshall  plasma 
armature  seem  to  indicate  that  conductivity  will  play  an  important  role  in 
determining  the  growth  rates  of  R-T  instabilities.  For  small  x,  the 
dimensionless  mass  density  is  of  order  x  1,  Examination  of  equations  (11) 
through  (14)  reveals  that  the  largest  terms  in  the  coefficients  A,  B,  C,  D,  and 
E,  are  of  order  x  * ,  making  them  several  orders  of  magnitude  larger  than  some 
of  the  smaller  terms  involving  a  a  dependence.  If,  because  of  this,  the 
coefficients  C  and  D  cannot  be  neglected,  the  dispersion  relation  will  depend 
on  the  solutions  of  a  fourth  order  ordinary  differential  equation  (ODE), 
instead  of  the  second  order  ODE  in  Tsai's  model. 

The  first  graph,  approximates  the  conditions  of  the  Westinghouse  EMACK 

rail-launch  device.  For  the  physical  parameters  reported  by  Deis  and  McNab**, 

again  employing  the  numerical  calculations  for  temperature,  electron  number 
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density,  ionization  state,  and  arc  length  of  Powell  and  Batteh  ,  the  EMACK 
rail-launch  device  can  be  characterized  by  x  =  0.004,  s  =  1.97  In  this  case, 
both  s  and  x  are  an  order  of  magnitude  smaller  than  for  the  Rashleigh-Marshall 
plasma  arc,  and  it  is  expected  that  finite  conductivity  effects  would  be  more 
pronounced. 

The  dispersion  graph  I  differs  markedly  from  the  ones  derived  under  the 
assumption  of  infinite  conductivity.  Growth  rates  of  short  wavelength 
perturbations  do  not  reach  arbitrarily  large  values,  but  appear  to  have  a 
finite  limit.  The  long  wavelength  behaviour,  however,  appears  to  be  similar 
irrespective  of  conductivity. 


While  the  growth  rate  of  the  R-T  instability  is  reduced  with  the  inclusion 
of  finite  conductivity  effects,  it  must  be  noted  that  even  the  longest 
wavelength  modes  permitted  by  the  geometry  of  the  gun  have  a  value  of  w  of 
order  1.0  (i.e.,  the  e-folding  time  of  the  disturbance  is  of  the  order  of  the 
time  it  takes  the  plasma  arc  to  accelerate  its  own  length).  Since  a  typical 
rail  device  is  about  30  times  longer  than  the  best  estimates  for  the  plasma  arc 
length,  there  is  ample  time  in  the  conrse  of  an  acceleration  for  R-T 
instabilities  to  set  in  and  develop  fully. 

A  data  point  was  picked  from  Graph  I,  and  the  associated  values  of  q  and  w 
were  inserted  in  equation  (10).  Numerical  integration  revealed  this  dispersion 
relation  to  be  for  an  n  =  1  mode  of  R.  Higher  order  modes  were  found  for 
various  valves  of  s  and  x,  but  they  had  smaller  growth  rates.  This  is  in 
agreement  with  the  work  of  both  Tsai  and  Powell,  who  documented  the  existence 
of  higher  order,  slower  growing  modes  under  the  assumption  of  infinite 
conductivity. 

Variations  in  s  and  x  noticeably  affect  the  growth  rates  derived  by 
numerical  integration  of  equation  (10).  Graph  II  gives  the  dispersion  relation 
of  the  n  =  1  mode  derived  for  the  values  s  =  0.1  and  x  =  0.15  on  the  same  scale 
as  Graph  I  for  comparison. 

Finite  conductivity  effects  are  even  more  pronounced  as  x  and  s  approach 
each  other  in  size.  Graph  II  shows  the  dispersion  relation  for  a  system 

characterized  by  x  =  1,  s  =  1.  While  the  growth  rate  diminished  for  shorter 

wavelength  disturbances,  it  is  interesting  to  note  that  the  real  part  of 
a)  grew  monotonically  to  an  asymptotic  value  of  about  4.3  (the  mode  found  for 
the  EMACK  device  was  non-oscillatory) . 

Integration  of  equation  (10)  for  a  point  from  Graph  II  reveals  that  this 
graph  is  also  for  an  n  =  1  mode.  Graph  III  plots  the  mass  current  density  R  as 
a  function  of  arc  position  for  this  mode. 

1  0 

Powell  and  Batteh  have  determined  how  temperature,  plasma  arc  length, 
acceleration,  etc.  in  a  rail  device  scale  with  current,  projectile  mass,  and 
other  factors.  Employing  these  scaling  laws,  and  assuming  a  Spitzer 
conductivity,  it  is  found  that  the  Rashleigh-Marshall  plasma  parameters  would 

scale  to  about  s  =  1,  x  =  1  for  an  initial  current  per  rail  height  of  roughly 


3  x  101*  statamp/cm.  While  this  may  not  be  feasible  at  present,  the  scaling 

relations  do  seem  to  indicate  that  the  growth  rate  for  R-T  instabilities  is 

diminished  slightly  with  higher  currents. 
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